Strong magnon softening in tetragonal FeCo compounds 
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Magnons play an important role in fast precessional magnetization reversal processes serving as 
a heat bath for dissipation of the Zeeman energy and thus being responsible for the relaxation of 
magnetization. Employing ab initio many-body perturbation theory we studied the magnon spectra 
of the tetragonal FeCo compounds considering three different experimental c/a ratios, c/a =1.13, 
1.18, and 1.24 corresponding to FeCo grown on Pd, Ir, and Rh, respectively. We find that for 
all three cases the short-wave-length magnons are strongly damped and tetragonal distortion gives 
rise to a significant magnon softening. The magnon stiffness constant D decreases almost by a 
factor of two from FeCo/Pd to FeCo/Rh. The combination of soft magnons together with the 
giant magnetic anisotropy energy suggests FeCo/Rh to be a promising material for perpendicular 
magnetic recording applications. 

PACS numbers: 71.15.Qe, 71.45.Gm, 75.30. Ds, 71.20.Be 
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Since the introduction of the first commercial hard disk 
drive in 1956, the recording density in a hard disk , that is 
the amount of information that can be stored per square 
inch, has increased by more than 7 orders of magnitude 
to meet an ever-growing needii This has been achieved 
by a simple scaling of the dimensions of the bits recorded 
in storage medium. Due to the superparamagnetic effect, 
however, the recording density has an upper limit. For 
longitudinal magnetic recording it is around 200 gigabit 
per square inch, whereas it is predicted to be much larger 
for perpendicular recording, up to 1000 gigabit per square 
inch, though this limit is constantly changing with the 
discovery of new materials j^^— 

The major problem in designing magnetic storage me- 
dia is to retain the magnetization of the medium over 
a long period of time despite thermal fluctuations. If 
the ratio of the thermal energy k^T to the magnetic 
energy per grain K^V, where V is the grain volume 
and K^ is the uniaxial magneto-crystalline anisotropy 
energy, becomes sufficiently large, the thermal fluctua- 
tions can reverse the magnetization in a region of the 
medium destroying the data stored therei^i^ In order to 
further increase the recording density in future record- 
ing media, high-i^u materials are neededji^ Additionally, 
a large saturation magnetization Mg is beneficial to re- 
duce the write field, which has to be applied by the writ- 
ing head. Materials that combine the desired large val- 
ues of Ku and Mg are tetragonal near-equiatomic FeCo 
alloys. The large values of K^ and Mg in these al- 
loys were first predicted by first-principles calculations^ 
and then confirmed by experiments!^^— In particular, 
Yildiz et ali^ achieved a strong perpendicular magnetic 
anisotropy (PMA) in tetragonal FeCo alloys epitaxially 
grown on Pd {c/a = 1.13), Ir {c/a = 1.18), and Rh 
{c/a — 1.24) substrates. The authors found that the 
PMA is very sensitive to the tetragonal distortion and 
increases with increasing c/a ratio, which allows to tune 
the PMA by growing the alloys on different substrates. 

Besides large K^ and Mg values, another very im- 



portant issue in magnetic recording applications is 
the magnetic switching time, which imposes physical 
limits on data rates and areal recording densitiesJ^ 
In current devices the switching speeds have reached 
a point where dynamical effects are becoming very 
importanti^"— Magnons are created in fast (field driven) 
as well as ultrafast (laser induced) magnetization rever- 
sal processes. ^^rJ^ The former case is of particular in- 
terest for current device applications. It is found that 
above some threshold magnetic field the uniform pre- 
cessional mode, i.e., the k = magnons decay into 
nonuniform magnons (k ^ 0), i.e., the Zeeman en- 
ergy stays in the magnetic subsystem and scatters be- 
tween magnon modesi^''— However, in ultrafast mag- 
netization reversal the high-energy electrons generated 
by the laser field decay into the lower-energy magnon 
excitationsi^""— In both cases spin-orbit coupling (SOC) 
is responsible for the transfer of the angular momen- 
tum to the lattice through different scattering mechanism 
like magnon-magnon, magnon-phonon, magnon-impurity 
scattering, and so on, where each process has a different 
relaxation time.^^'^^ The Landau-Lifshitz-Gilbert (LLC) 
equation with a phenomenological damping constant a is 
commonly employed to describe magnetization dynamics 
of small-angle precessional switching. ^^ However, recent 
studies have shown that in the case of large-angle (fast) 
switching, in which the magnons are created, the LLC 
equation should be extended in several aspects^i^, in 
particular a k-dependent damping constant a^ has been 
proposed )22ii2^ which allows short-wavc-lcngth magnons 
to relax faster than those with k — > 0. Thus, the magne- 
tization relaxation processes, specifically the damping of 
magnons, play an important role in designing ultrahigh- 
density magnetic recording media. 

The aim of the present Communication is to study 
magnon dynamics in tetragonal FeCo compounds from 
first principles. Using a recently developed Green- 
function method'^'* based on the many-body perturbation 
theory in the GW approximation in combination with the 



TABLE I: Lattice parameters a, spin magnetic moments m° (in fis), average screened on-site direct (diagonal) {W = 



1 ^(^3d) \Y„n;nn) and exchange {J ^ ^Yl 
magnon stiffness constants D (in meVA^ 
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) Coulomb matrix elements (in eV) between the 3d orbitals, and 
for tetragonal FeCo compounds grown on Pd, Ir, and Rh. Lattice parameters are 



taken from Ref.[T3. 
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multiple-scattering T matrix in a Wannier basiS) i we 
have calculated the dynamical spin susceptibility (DSS) 
of tetragonal FeCo compounds whose c/a ratios were 
fixed to the experimentally determined values that relate 
to the three different substrates. As the unit cell con- 
tains two magnetic atoms, the calculated magnon dis- 
persions exhibit two branches: an acoustic and an op- 
tical branch. The former persists throughout the Bril- 
louin zone indicating a localized nature of magnetism in 
FeCo compounds. The optical branch, on the other hand, 
is heavily damped due to the coupling to single-particle 
Stoner excitations. We find that the tetragonal distortion 
gives rise to significant magnon softening. The average 
magnon stiffness constant D decreases almost by a fac- 
tor of two from FeCo/Pd to FeCo/Rh, which means that 
acoustic magnons can be excited much more easily in the 
latter material than in the former one. In field-driven fast 
magnetic switching processes, which take place on a time 
scale of ns to 100 ps, the excess Zeeman energy will be 
transferred to the acoustic magnons and thus magnon 
stiffness constant D and life time of k 7^ magnons play 
a decisive role determining the strength of the switch- 
ing field and switching time. The latter is limited by 
the damping of magnons. Furthermore, damping also 
prevents "back-switching" effect, which reduces the data 
rates in magnetic recording devices4i 

To calculate the ground-state properties of the tetrag- 
onal FeCo compounds we use the full-potential linearized 
augmented plane-wave method as implemented in the 
FLEUR code<^ The exchange-correlation potential is cho- 
sen in the generalized gradient approximation.'^^ The 
muffin-tin radii of the Fe and Co are chosen to be 2.29 
a.u. A dense 16 x 16 x 16 k-point grid is used. The max- 
imally localized Wannier functions are constructed with 
the Wannier90 codei^i^ The DSS is calculated within a 
T-matrix approacb^^ as implemented in the SPEX codeSi 
using 8000 k-points in the full Brillouin zone. We briefly 
review the method here. Within many-body perturba- 
tion theory the transverse DSS, x ^j can be schemati- 
cally written as x ^ = Xks^ + Xks -^ ^Xks ' where the 
first term on the right-hand side represents the response 
of the noninteracting system, i.e., the Kohn-Sham DSS. 
The second term contains the T matrix, which is given 
by T '" = [1 — W^Xks"]""^^' where W is the screened 
Coulomb interaction. The T matrix describes dynami- 
cal correlation in the form of repeated scattering events 
of particle-hole pairs with opposite spins and is respon- 



sible for the formation of collective magnon excitations. 
Details of the formalism, implementation, and applica- 
tions to 3d transition metals can be found in Ref. 35|. 
The DSS provides complete information on the magnetic 
excitation spectrum including collective magnon modes 
as well as single-particle Stoner excitations together with 
their respective lifetimes4ir— We note that magnon life- 
times and Stoner excitations are not accessible within the 
adiabatic approximation, a method mostly employed so 
far for the calculation of the magnon dispersion within 
density functional theory.^S 

Experimentally, FeCo alloys have been grown on the 
Pd, Ir, and Rh substrates in the body-centered tetrag- 
onal structure, in which the in-plane lattice constant is 
determined by the substrate and the out-of-plane lattice 
constant changes so as to keep the volume constantJ^ 
The experimental lattice parameters used in the calcula- 
tions are presented in TablelH Yildiz et al. have shown 
that the structure remains stable for film thicknesses of 
up to 15 monolayers, which is large enough to consider it 
as a bulk in the context of theoretical modelling.''^ Since 
no experimental information is available on the micro- 
scopic atomic order of FeCo alloys grown on the different 
substrates, we assume a tetragonally distorted CsCl-type 
(B2) structure derived from the known cubic bulk phase. 
We note that the mechanism behind the giant uniaxial 
MAE observed in tetragojial FeCo compounds has been 
discussed in detail in Ref.|7 and will not be dwelt on here. 
Indeed, our calculated values of uniaxial MAE (results 
not shown) are very similar to those reported by Burkert 
et al7 

We start with a discussion of the magnetic moments 
and the matrix elements of the screened Coulomb po- 
tential W . The latter are a crucial ingredient for the 
construction of the DSS. The calculated values for the 
three different c/a ratios are presented in TableU) As 
seen, the spin magnetic moment of the Fe sublattice is 
substantially enhanced with respect to bulk bcc (or fee) 
Fe, which has a magnetic moment of about 2.2 /i^, while 
the Co sublattice shows moments that are more simi- 
lar to the corresponding value of bulk Co, 1.62 jib- In 
contrast to bulk Co, bcc Fe is a weak ferromagnet, and 
thus its magnetic moment is very sensitive to the local 
environment. The total spin magnetic moment of the 
unit cell is around 4.4 /i^ and almost insensitive to the 
tetragonal distortion, which stems from the strong fer- 
romagnetic nature of the FeCo compounds [see inset in 
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FIG. 1: (Color online) (a) Imaginary part of the un- 
renormalized Kohn-Sham spin susceptibility of the tetragonal 
FeCo compounds grown on Pd {c/a = 1.13) for selected wave 
vectors. The inset shows the spin-resolved total density of 
states; (b) the same for the renormalized spin susceptibility. 
Note the different scales of the axes in the two figures. 



Fig-UJa)]. Such a large magnetic moment is desirable for 
magnetic recording applications as it reduces the write 
field of the writing head. Not only magnetic moments 
but also the average screened on-site Coulomb matrix el- 
ements W (direct) and J (exchange) of FeCo compounds 
are insensitive to the tetragonal distortion. The obtained 
values are slightly larger than the corresponding values 
in the bulk phase of bcc Fe and fee Co. As discussed in 
detail in Ref.|35|, this difference can be attributed to the 
larger exchange splitting of the Fe and Co atoms in FeCo 
compounds as the larger the exchange splitting the less 
screening takes place leading to a stronger Coulomb inter- 
action W . Indeed, with increasing c/a ratio the magnetic 
moment (exchange splitting) of the Co atom decreases 
slightly, giving rise to a small reduction in the Coulomb 
matrix elements Wco as shown in TablelH 

Figure [TJ^a) presents the non-interacting Kohn-Sham 
transverse DSS of the tetragonal FeCo compounds grown 
on Pd [c/a = 1.13) for selected wave vectors at high sym- 
metry points M [q = (0.5, 0.5, 0)], Z [q = (0, 0, 0.5)], and 
R [q = (0.5, 0, 0.5)]. As there is no dynamical correlation 
due to the absence of electron-electron interactions, only 
single-particle spin- flip Stoner excitations exist. As a 
consequence, the spectral function Imx^g (q, w) exhibits 
a broad peak originating from spin-flip transitions be- 
tween occupied majority and unoccupied minority states. 
The peak maximum at around 2.5 eV corresponds to the 
exchange splitting of the FeCo compounds visible in the 
density of states shown in the inset of Fig.[lja). The sit- 
uation is very similar for FeCo compounds grown on Ir 
and Rh. As will be discussed below, Stoner excitations 
are responsible for the damping of the magnons. 

When dynamical correlation is included in the form of 
interacting spin-up and -down electrons via the screened 
Coulomb interaction W (see TableHl , additional magnon 
peaks appear in the spectral function of the interact- 
ing system Imx '"(q, cj) at low energies as illustrated in 
Fig-IHb) for the case of FeCo/Pd. Since the unit cell con- 
tains two magnetic atoms, we obtain two modes: a low- 
energy acoustic mode and a high-energy optical mode. 
The former corresponds to the in-phase precession of the 
Fe and Co magnetic moments, while the latter is due to 




r Z R A 



FIG. 2: (Color online) Magnon dispersion of the tetragonal 
ordered bulk FeCo alloy as grown on (a) Pd {c/a — 1.13) (b) 
Ir {c/a = 1.18), (c) Rh {c/a — 1.24), along the high-symmetry 
lines in the Brillouin zone. 



out-of-phase precession. The broadening of the peaks is 
caused by coupling to single-particle Stoner excitations. 
Plotting the renormalized susceptibility as a function of 
the wave vector yields the magnon dispersion as displayed 
in Fig. [2] along the high-symmetry lines in the Brillouin 
zone for all three systems. As seen in all cases, we obtain 
a well-defined low-energy acoustic branch throughout the 
Brillouin zone, while the high-energy optical branch lies 
above 0.75 eV and is heavily damped due to coupling to 
Stoner excitations. This strong damping can be traced 
back to the density of states (DOS) of Stoner excitations 



in the respective energy region, see Fig.[lja) . We observe 
a drastic increase of the DOS above 0.8 eV. Below 0.8 
eV the DOS is small because the FeCo alloy is a strong 
ferromagnet, i.e., the spin-majority states are fully occu- 
pied. Thus, though damped, the optical branch remains 
identifiable in FeCo alloys, while in weak ferromagnets 
like bcc Fe the magnons persist only up to 500 meV. 
Above this energy they disappear due to strong coupling 
to Stoner excitations4^i^ We note in passing that s-d 
mixing leads to damping of the magnons also in cases of 
strong ferromagnetism, but this effect is relatively weak 
compared to damping in weak ferromagnets. Damping of 
magnons does not mean that the angular momentum is 
transferred to the lattice. In the absence of SOC the an- 
gular momentum stays in the magnetic subsystem, i.e., it 
is transferred from magnons to the non-coherent single- 
particle spin-flip Stoner excitations. The SOC is respon- 
sible for angular momentum transfer from the magnetic 
subsystem to the lattice subsystem. We also note that 
we have not included SOC in the calculation of the DSS. 
The SOC gives rise to an opening of a gap in the magnon 
dispersion at the F point, whose value is determined by 
the MAE Ka, which is less than 1 meV for the systems 
considered here.— Apart from that, we do not expect a 
qualitative change in the magnon spectra because the 
SOC only has a negligible effect on the electronic struc- 
ture. 

With increasing c/a ratio the two branches are pushed 
in opposite directions, i.e, the excitation energy of the 
optical (acoustic) magnons increases (decreases). This 
behavior, on the one hand, leads to weak damping of the 
acoustic magnons since the intensity of the Stoner DOS 
decreases at low energies, while the situation is just the 
opposite for the optical magnons. On the other hand, 
it gives rise to magnon softening, i.e., the magnon stiff- 
ness constant D presented in TableU is considerably re- 
duced with increasing c/a ratio. As seen, from FeCo/Pd 
to FeCo/Rh the average D decreases almost by a fac- 
tor of two. The calculated average magnon stiffness 
constants D for FeCo/Pd and FeCo/Rh are close to 
the experimental values of fee Co {D = 580 meVA^) 
and bcc Fe {D — 280 meVA^), respectively!^ Further- 
more, the in-plane (D|| ) and out-of-plane (D^) magnon 
stiffness constants differ a lot, and this difference in- 
creases with increasing c/a ratio. The small value and 
anisotropy of the exchange stiffness constant D can be 
attributed to a strong direction dependence of the ex- 
change interactions.-- Note that in itinerant ferromag- 
nets there are several coexisting exchange interactions. 
A detailed discussion of them for present systems is be- 
yond the scope of the this work. In the following, we will 
discuss them only qualitatively. In 3d ferromagnets and 
their alloys the total exchange coupling can be divided 
into two contributions: Jt = ^direct + -/indirect, where 
the first term (direct coupling) is a short-range interac- 
tion due to the overlap of the 3d wavefunctions and its 



strength depends on the distance between the magnetic 
atoms, while the long-range indirect part is due to the 
coupling of the localized 3d moments to the itinerant sp 
electrons. For c/a ^ 1 the Jdircct becomes anisotropic. 
With increasing c/a, i.e., from FeCo/Pd to FeCo/Rh, 
the in-plane and out-of-plane components of Jdircct can 
differ a lot. The former (latter) is expected to increase 
(decrease) due to smaller (larger) interatomic distances. 
Consequently, the c/a behavior of the direct exchange 
coupling can qualitatively account for the anisotropy of 
the magnetization and the reduction of the magnon en- 
ergies (acoustic branch) along the z direction in FeCo 
compounds. However, the strong in-plane magnon soft- 
ening is more likely connected to c/a behavior of the long 
range indirect exchange interactions, which give a sub- 
stantial contribution to the total exchange coupling Jt 
with a negative sign. A qualitative estimate of its contri- 
bution to Jt is not easy without very detailed electronic 
structure analysis since this coupling shows RKKY-type 
oscillations, extends over very large distances, and is very 
sensitive to tetragonal distortioui^ Its strength and long- 
range behavior is determined by several parameters such 
as conduction electron spin polarization, Fermi surface 
topology, position of unoccupied states with respect to 
the Fermi level and so on. For a detailed discussion on 
the indirect exchange coupling in 3d transition metal al- 
loys the reader is referred to Ref.i51j. Finally we would 
like to note that as the magnetism in itinerant ferro- 
magnets depends on the electronic states far from the 
Fermi level, the disorder between Fe and Co sublattices is 
not expected to substantially influence the magnon spec- 
tra of FeCo compoundsi-^-^ However this is not the case 
for MAE, which is very sensitive to the Fermi surface 
topology.^'^ 

In conclusion, we have calculated the magnon spectra 
of the tetragonal bulk FeCo compounds from first princi- 
ples considering three different experimental c/a ratios: 
FeCo grown on Pd, Ir, and Rh with c/a =1.13, 1.18, and 
1.24, respectively. We have found that for all three cases 
the short-wave-length magnons are strongly damped and 
tetragonal distortion gives rise to a significant magnon 
softening. The magnon stiffness constant D decreases 
almost by a factor of two from FeCo/Pd to FeCo/Rh, 
which reduces the switching field and yields efficient exci- 
tation of the k 7^ magnons. Furthermore, the obtained 
strong damping of large-wave-vector magnons in FeCo 
compounds suggests a k-dependent damping constant 
ak in the LLC equation in describing magnetization dy- 
namics of large-angle fast precessional switching. Com- 
bination of soft magnons with their substantial damp- 
ing at large wave vectors as well as giant MAE suggests 
FeCo/Rh to be a very promising material for ultrahigh- 
density perpendicular magnetic recording applications. 
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